Aspects of Graviton Detection: Graviton Emission 
and Absorption by Atomic Hydrogen 

Stephen Boughn* and Tony Rothman^ 

* Haverford College , 
Haverford PA, 09140 
f Princeton University, 
Princeton, NJ 08544 

MfeX-ed February 6, 2008 



Abstract 

Graviton absorption cross sections and emission rates for hydrogen are cal- 
culated by both semi-classical and field theoretic methods. We point out several 
mistakes in the literature concerning spontaneous emission of gravitons and re- 
lated phenomena, some of which are due to a subtle issue concerning gauge 
invariance of the linearized interaction Hamiltonian. 
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1 Introduction 



Two years ago, Dyson published a conjecture that no conceivable experiment per- 
formed in our universe could detect a single graviton pQ. Recently, in a companion 
paper [2] (henceforth RB), we have addressed Dyson's proposition and have been un- 
able to find any clear-cut fundamental reason forbidding the detection of one graviton. 
However, when anything remotely resembling "real" physics is taken into account, a 
detection becomes impossible, making Dyson's conjecture very likely true, at least 
without the introduction of exotic physics, such as extra dimensions. 

*sboughn@haverford.edu 
t trothman@princeton.edu 
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In the process of checking our calculations against graviton spontaneous emission 
rates for hydrogen, we have discovered several mistakes in the literature involving both 
spontaneous emission rates and the response of elastic media to classical gravitational 
waves. Although some of the mistakes are evidently numerical, others involve a rather 
subtle issue of gauge invariance of the interaction Hamiltonian in linearized gravity. 
Our purpose here is to present these calculations in detail, thereby illuminating the 
source of the difficulties. To demonstrate the validity of our results, we compute 
transition rates from the 3d to Is hydrogen states using both semi-classical and field 
theoretic methods in two different gauges. We also present the consistency check 
that led to the discovery of the mistakes: the requirement that detailed balance be 
satisfied among the spontaneous emission, stimulated emission and absorption rates. 
Finally, we compute the graviton ionization cross section of hydrogen, because in RB 
we chose this gravitational analogy to the ordinary photoelectric effect as the most 
relevant one in addressing Dyson's question. 



2 Semi-classical Analysis of Spontaneous Graviton 
Emission 

In the standard semi-classical treatment of spontaneous emission of photons from 
atomic dipole transitions (see, for example, Schiff the transition rate is found by 
beginning with the classical expression for power emitted by a dipole, 

P=^-f\M\ (2.1) 

3c 

in which J = / 3(r)d 3 r is the total current and k is the wave number of the emitted 
radiation. One then attempts to interpret this expression in a quantum mechan- 
ical way by identifying the current density J (charge density times velocity) with 
the probability current e|\l/| 2 p/m, for wavefuction \1/ and electron charge e. Be- 
cause the quantum mechanical transition is between states two states a and b, the 
standard treatment replaces e\^\ with e^*^, and so the total current becomes 
Jo = — / ty*p^i,d 3 r. Finally, the quantum mechanical transition rate is assumed to 
be 

r _ P 

with P given by inserting the "quantum" expression for Jo into Eq. (|2.1jl . The 
justification for this argument is, to say the least, tenuous; therefore it is important 
to verify the resulting transition rate in some other way, e.g., that it satisfies detailed 
balance or that a "proper" field theoretic calculation yields the same result jSj- 

For the case of spontaneous graviton emission, we begin with the (quadrupole) 
gravitational analog of Eq. ([2.1)1 . 

p = -mtUM- ( 2 - 2 ) 
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In this expression, Ij k = J p(r)(xjXk — ^5jkr 2 )d 3 r is the reduced mass quadrupole 
moment, r = {xj} are cartesian coordinates in the local inertial frame, and p is the 
mass density of the radiating system J3J. The normalization for Ijk, as well as J above, 
is usually taken to be Ijk(r,t) = Ijk(r)e~ tuJt + c.c. Following the electromagnetic 
procedure, the semi-classical approach for gravity substitutes the electron mass m e for 
e and replaces p by m c $^ a . (The equivalent mass density associated with the electric 
fields of the electron and proton is, for non-relativistic systems, small compared to 
the electron mass; we discuss this issue in more detail in £0) Therefore the reduced 
quadrupole moment becomes 

Ijk = m eJ %(xjX k - ^6 jk r 2 )y a d 3 r. (2.3) 

To compare with previous work, we can now compute the transition rate from the 3d2 
to the Is state of hydrogen. In terms of the Bohr radius a = h /m e e 2 , the normalized 
wavef unctions are 

Inserting these expressions into into Eq. ()2.3|) and subsequently into Eq. ()2.2|) yields 
a transition rate of 

= p/nu = ^Grnl-y = cPGmfc = ^ x 

1 5 x 2 13 hc 5 360h 2 v ' 

where a = e 2 /he is the fine structure constant. It is straightforward to show that the 
transition rates for all 3d states (3c?(±2), 3d(±l), and 3<i0) are all the same. 

In his standard text Gravitation and Cosmology [3], Weinberg presented a calcu- 
lation of this transition rate using the same semi-classical approach. His result was 
T = 2.5 x lCT 44 ,^ 1 , which differs from ours by more than four orders of magnitude. 
Since the methods are identical, and the dimensional factors agree, the discrepancy 
is presumably due to a numerical error. 

Indeed, in his book on quantum gravity jSj, Kiefer obtained the result in Eq. 
(|2.5|) . again using a semiclassical analysis. While this agreement may indicate the 
above analysis has no computational errors, the semi-classical treatment itself "can 
claim only a moderate amount of plausibility" 0- Stronger support is afforded by 
demonstrating that the transition rate of Eq. (|2.5|) satisfies detailed balance. To do 
this it is first necessary to evaluate the cross sections for absorption and stimulated 
emission. 



3 Absorption, Stimulated Emission, and Detailed 
Balance 

Unlike spontaneous emission, the processes of absorption and stimulated emission 
can be treated adequately in the context of classical quantum mechanics, that is, 
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without the use of quantum field theory. According to first order perturbation theory, 
the transition probability between two hydrogenic states ty a and is proportional 
to the square of the matrix element (^/b\H\^/ a ), where the interaction Hamiltonian 
H is derived from the interaction Lagrangian by its definition H = pv — L. For 
metric deviation h^ v = g^ v — rj^ v « 1 and stress-energy tensor T^, the interaction 
Lagrangian density is given by Ej 

£ = hi, u T^. (3.1) 

Although the action from which £ is derived is a scalar, £ itself is not gauge invariant. 
It is in fact this property that has led to the more serious conceptual errors mentioned 
in the Introduction and which we will discuss in detail in £0 The standard choice 
is to work in a local inertial frame (LIF) in which case the dominant term of the 
stress-energy tensor is just the mass-energy density, and so £ ~ ^h 00 T 00 . When the 
generalized velocities are negligible, as in an LIF, the interaction Hamiltonian density 
is simply 7i = —£, and the corresponding interaction Hamiltonian for a localized mass 
m is therefore H m — |m/ioo> if we assume that /ioo( r ) is approximately constant in the 
region of nonvanishing mass density. Although LIFs are the most nearly Minkowskian 
and hence physical frame of reference, it is in the transverse-traceless, or TT, gauge 
that gravitational waves are most easily interpreted. 

The TT gauge is a subset of what is variously called the Hilbert, Einstein, Fock, 
de Donder or Harmonic gauge, which is analogous to the Lorentz gauge of electromag- 
netism and defined by the requirement that h 1 * Utfi = The TT gauge makes the 

further choice that h^ = ^ = 0. Then for an amplitude h and polarization tensor 
e M „ the metric deviation can be written as h^ u = he^ v , where e o = e^o = M = 0. 
Consequently, a harmonic, plane gravitational wave (GW) can be expressed as 

h% = he* k "*-« e >e Jk + c.c. (3.2) 

where j, k, I = 1, 2, 3 and u = kc. In this case the Hilbert condition becomes ekik 1 = 0, 
which for a plane wave propagating in the z-direction gives e xx = —e yy = e xy = e yx 
and e Z j = 0. In what follows, we normalize the nonzero components of the polarization 
tensor to |ey| = l/\/2. 

Independent of normalization, /i o i n a LIF can be expressed in terms of the 
quantities hj£ of the TT gauge. Since the duration of the interaction of the incident 
GW with the physical system (the hydrogen atom in this case) may be long, it is 
necessary to use a coordinate system that is locally intertial for an extended time. 
There are an infinite number of such coordinate systems; however, it can be shown 
that to 2nd order in coordinate displacement h 0Q is the same in all such systems. 
Perhaps the most useful of these are "Fermi normal coordinates" in which j3] 



h 00 — —R j f : x : 'x k 



(3.3) 
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where R a /3^s is the Riemann curvature tensor. To lowest order in h a p 

Ra(3~fS — -{h a S,p~( + ^/3 7 ,a5 — hp6,a-y ~ h al> ps) (3-4) 

and is gauge invariant. Combining Eqs. (JS2J), ()3.3|) . and ()3.4j) yields 

/i 00 = -~u 2 he i{ - kixl -^x j x k e jk + c.c. 
2 

The LIF interaction Hamiltonian therefore becomes 

H = -m e u 2 hx j x k e jk e^ hxl -^ + c.c. (3.5) 
This is the Hamiltonian we will use in computing the matrix elements. 

To calculate the transition rate, one also needs the flux of incident GWs, which is 
given by jl] 

where () indicates an average over several cycles. Assuming equal amplitudes for the 
two polarizations and the normalization |e^| = 1/Vz, we have 

T = — -h 2 . 

8nG 

The transition probability per unit time between two discrete states is not constant 
in time if the incident radiation field and Hamiltonian are strictly monochromatic. So 
we assume that the radiation is spread over a range of frequencies with uncorrelated 
phases. In a small neighborhood Au of each of these frequencies the flux is given by 
dT = I(u)Au where is the GW intensity. Therefore, 

h2= 8kG ^ 
The transition probability per unit time is then given by standard perturbation theory 

m 

r 1 ^ 4|(fr m |ij-|fr n )| 2 sin 2 \(oj mn - u)t 

t ~ h 2 (uj mn - id) 2 

_ 27iGm 2 e r ; u 2 I(u)Au sin 2 \{uj mn - u)t 



(3.7) 



Here, the difference between the two electron energy levels is AE = Two mn and the 
sum is over the frequencies of the phase independent GWs each within a band Au>. 
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In the limit that Aoo is infinitesimally small, the summation can be replaced by an 
integral. Since the time factor has a sharp maximum at u = u> mn , I{uj)/uj 2 can be 
taken outside the integral and the limits extended to ±00. Eq. ()3.7|) then reduces to 

r = ^ G <^nJ{u mn ) | J ^^^^.^^2 

Finally, for the transition from the 3d to Is states of hydrogen, the wavelength of 
the GW is much larger than the extent of the wavefunction so that e tklX ~ 1 (the 
"dipole" approximation). Then 



Tx 2 Gm 2 e u? mn I{b) mn ) 



J ^* m x j x k e jk ty n d 3 r\ 2 



where 



^2^3 \- L 'jk e jk\ > \ 6 -°) 



Dij = J ^ b XiXj^ a d 3 r. (3.9) 



For a GW propagating in the z direction, the two linear polarizations are as given 
above and 

\D jk e jk \ 2 = 2\J ^ m xy^ n d 3 r\ 2 + i| J ^ m (x 2 - y 2 )V n d 3 x\ 2 . (3.10) 



To compute the mean transition rate for GW's incident from all directions we 
simply average \Dj k ej k \ 2 over the sphere. The following is due to Dyson [§]. Since 
electron wave functions do not depend on the direction of the gravitational field, we 
can introduce two orthogonal unit vectors A, ft and rewrite Eq. (j3.1()|) as 

\D jk e jk \ 2 = 2\J ^* m \ ■ r/i • r^ n d 3 r\ 2 + ^\ J V* m (X ■ r A • r - jl ■ rjl ■ r^ n d 3 x\ 2 . 
The average of this expression over all directions is then 

< \Dj k ej k \ 2 > = ^ ^ dQ,\^\ijljDij\ k (iiD* kl + Xi\jDij\ k \iD* kl 

1— . . 1.. — ( 311 ) 

— ^XiXjDijfikfiiD^ — jfL k fLiD k i\iXjD*j], 

where repeated indices are summed and 

< XiXjDijX k XiD* kl > = < fiifijDijfikfiiD^ > . 
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Because p, and A are orthogonal, we can eliminate \x by the following trick: Pick a 
direction for A, say A = k. Then in the usual spherical coordinates ft = icos<f)+jsin<j). 
The average < fi k fii > over the unit circle in the plane perpendicular to A can be seen 
to be < fi k fii >= 1/25^ — l/2AfcA;. Since this is a tensor equation it is true in any 
coordinate system and, hence, for any direction A. Inserting this expression into Eq. 
(J3.11j) and making use of the identity 

. . . . If .... 1 

< A^AjAfcA; >= — J dft AjAjAfcA; = — (^ij^kl + ^ik^jl + fiilfijk) 

yields 

< \D jk e jk \ 2 >= \{D ]k D* k - l -D 03 Dl k ). (3.12) 

It is this average that must be used in Eq. ()3.8j) to compute V. For the transition 
from the Is to the 3d2 state of hydrogen, the appropriate values of of Dj k are, from 
Eqs. (El and Q, 



D Z z D xz Dy Z 



4„2 



3 4 a 

D X x ^yy iD X y ^ . 

Substituting these values into Eq.f I3.12"j) yields 

8 3 8 a 4 
<\D ]k e jk \ 2 >=-\D xx \ 2 = ^ ¥ - r (3.13) 

Finally, the absorption rate for the Is to 3d2 transition follows from Eq. ()3.8|) : 

r _ S 8 7i 2 Gm 2 e a 4 ul n I(u mn ) 

5 x 2i3^ 2 c 3 ^' i4J 

This expression also gives the stimulated emission rate between the 3d2 and Is hy- 
drogenic states. 



To check the consistency of the absorption rate with the spontaneous emission rate 
calculated in ^21 we assume that gravitational radiation within a cavity is in thermal 
equilibrium with emitters and absorbers in the walls at temperature T. (While this 
detailed balance argument is "valid" in a certain sense, it is straightforward to show 
that no such cavity can in principle be constructed for gravitational waves.) Equating 
the radiation absorbed per unit time with that emitted per unit time we have 

N3d2(X sp + T st ) = N ls T ab , 

where T sp , T st , and T a b are the transition rates for spontaneous emission, stimulated 
emission, and absorption, is the number of atoms in the 3d2 state, and iV ls is the 
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number atoms in the Is state. Following Einstein, we expect N^/Nu = e~ nu;mn ^ kT . 
Now, substitute Eq. (|3.14|) and the spontaneous emission rate from Eq. ()2.5|) into the 
above expression. Solving for I(u mn ) yields 

,3 

which is consistent with the intensity of black body radiation. This indicates that 
the ratio T sp /T ab is correct and that the result yielded by the semi-classical approach 
in glis as valid as Eq. (jTTIjl . 

The absorption rate Eq. 1)3.14)1 can be expressed in terms of an integrated cross 
section, 

1 



r = — / <r(u)I(u)doj. (3.15) 

Clearly the absorption cross section is sharply peaked near u mn , and so / a(uj)duj = 
KojT J I(u mn ). If we define an average cross section as < a >= J a(u mn )duj /ui mn , then 
from Eq. (jSHj): 

3 8 n 2 GmWuL, 

< o >= 



5 x 2 13 hc 3 

For the transition between the Is and 3d2 state u mn = yielding 

3 4 vr 2 Gh „ 2 . 

aahs = = °- 3Upl - (3 - 16) 



Cr 



Here, £ p i ~ 10~ 33 cm is the Planck length. Surprisingly, all the physical constants 
associated with the hydrogen atom have disappeared from the cross section, leaving 
only the square of the Planck length and a numerical constant of order unity. We 
return to this important point in §H1 

The absorption rate (j3.14j) was for unpolarized GWs averaged over all incoming di- 
rections. For completeness, the absorption rate for one polarization of a gravitational 
wave incident in the 9, direction is readily shown to be 

3 8 7r 2 Gmla 4 u!l in I(uj mn ) r 2 2 2 2 . 2 

T = — — q — 1 + cos 0) cos 26 + Acos Osm 2 

2 16 h 2 c 3 1 

The rate for the other polarization is obtained by interchanging the sin 2 2(j) and cos 2 2cf). 

4 Field Theoretic Calculation of Spontaneous Gravi- 
ton Emission 



As an independent check of Eq. (J2.5)) we next compute the transition rate via a 
field theoretic approach, that is, in terms of gravitons. This also allows a comparison 
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with the result of Lightman et al., who in problem 18.18 of their well-known Problem 
Book in General Relativity and Gravitation [TU| also used field theoretic methods to 
compute the 3d0 to Is transition rate in hydrogen. 

To quantize gravitational waves in the linearized theory we follow the standard 
procedure of decomposing the metric perturbations into plane waves: 

hjk = 4r E K^e^-^ + c.c. (4.1) 
v V M 

Here, h^ a are the Fourier amplitudes, a indicates the polarization, is the polar- 
ization tensor, and box normalization with volume V is assumed. The energy density 
of GWs corresponding to the flux T given by Eq. ()3.6|) is 



c 



2 



dE/dV = —{h^h 1 ^). 

Substituting Eq. ()4.1|) into this expression and integrating over all space gives the 
total energy in the GWs. Noting that 



and 



then 



( J e i(k.r-ort) e -i(k'.r-a/t) d 3 r j = V 5^5 U 



2 2 

& — \ n n e jk e jk ft k,a^k )Q ,'- 



16nG 



For the normalization |ejy| = l/\/2, we have J2jk e jt e jk" = &ota'i and 

c 2 oj 2 

E =l^G~P hkj2 - 

Next identify h^ a and h^a with the raising and lowering operators, a[_ a and ak, a , 
such that the classical energy and quantum Hamiltonian agree with one another, i.e., 

E = H = \ S huj ( a i,a a ^a + ak, a a ] Ka ) = 2~2( N ^« + ^) huJ 

k,o k,a — 

where N^ a = a a^ a is the number operator. One can always reset the zero of the 
energy scale so that H = £} k a ^-^k,a- The only way that E = H is if we make the 
identifications 



1 IQirGh 



c V oj 
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and 



Then 



1 167tGh t 



V* -> -777 \ —— W^e)C exp i(k • r - ut) + 4 e^ Q exp -z(k • r - ut) 
cVV k .„ V UJ L 



(4.2) 



We now consider the transition rate between two hydrogen states that involves 
the emission of a single graviton. According to the golden rule the transition rate per 
solid angle is given by 

§ = 2 i\W(f\H\iM\ 2 P- («) 

Here, \i)\0) is the initial electron and graviton state (with no gravitons), is the 

final electron and graviton state (with 1 graviton), and 

Vuj 2 

P 



(2n) 3 hc 3 

is the graviton energy density of states per solid angle, again box normalized. From 
Eqs. (|3.2|) and (|3.5|) . the interaction Hamiltonian (in the local inertial frame) can be 
written 

H = T I^h jk x j x k (4.4) 

where hjk is expressed in the TT gauge. In order to compare directly with the results 
of Lightman et al.([TU]) we compute the spontaneous transition between the 3d0 and 
Is states of hydrogen, although as pointed out earlier, the transition rates from all 
the 3dm states are identical. Substituting Eq. (|4.2j) into Eq. (|4.4|) the transition 
matrix element becomes 



(l|(ls|tf|3d0}|0> = ^f-^Jd 3 r% d0 e^x^ ls (4.5) 
where uj = (E 3 do — Eu)/h. Substituting Eq. (|4.5jl into Eq. (|4.3j) yields 



dQ iirhc 5 



where is defined in Eq. (|3.9j) . As before, \Djkejk\ 2 should be averaged over all 
directions. With Eq. (j3.12|) we then have: 

2m 2 e u 5 G 1 . 

T -^hc^ {D ^ k ~3 D " Dkk) - 
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For the 3d0 and Is hydrogen wavef unctions we find that D xx = D yy = —\D 
— a 2 w3 7 /2 15 and, therefore, 



r 38Gm ' flV (46) 

which is precisely the same as the semi-classical result in Eq. (|2.5|) . It is this agree- 
ment that provides the strongest confirmation of the semi-classical treatment. 

The Problem Book authors obtained a result that is larger than this one by about 
an order of magnitude. In part the disagreement is due to numerical and normaliza- 
tion errors. Nevertheless, when these mistakes are corrected their result still differs 
from Eq. (|4.6jl . The underlying reason can be traced to the fact that whereas our 
calculation was carried out in the local inertial frame, Lightman et al.pU] worked 
entirely within the TT gauge. One can recover the proper result in the TT gauge; 
however, much care must be taken to correctly interpret the interaction Hamiltonian 
in that gauge, as will become apparent in the following section. 



5 Gauge Properties of the Interaction Hamilto- 
nian 

In his pioneering work on gravitational wave detection Weber JT] used the equation 
of geodesic deviation to deduce the gravitational force density of weak gravitational 
fields acting on non-relativistic matter. In a local inertial frame the gravitational 
force density /J is 

ft = - P Rl ko x k , (5.1) 

where p is the mass density of the detector, and R 3 0k0 are components of the Riemann 
curvature tensor. Then the equation of motion of the detector is pd 2 x^ /dt 2 = /•? + f 3 g 
where fi is the total classical (i.e., non-relativistic) force density on the detector mass 
elements. This relation implicitly assumes that the gravitational wave has negligible 
effect on the physics of the detector, which is a reasonable assumption because the 
LIF is the most nearly Minkowskian frame. Since gravity couples to the mass-energy 
density of matter, one expects that the gravitational modification of the classical 
forces in the detector will be proportional to the binding energy density of the system. 
In fact, it can be shown that this is the case jS] . While Eq. (J5.1|) strictly holds true only 
in a LIF, i?Q fc0 is already first order in h^, and so the right hand side of the equation 
is, to first order, invariant under infinitesimal coordinate, or gauge, transformations. 
It is in this sense that the expression for the gravitational force in Eq. (j5.1|) is gauge 
invariant. 

As an alternative formulation to Weber's, Dyson [7j introduced the interaction La- 
grangian density, C = hh^T^" , which we used in ^21 and 2] above. It follows directly 
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from the general definition of the stress-energy tensor as the functional derivative of 
the matter action, J m , with respect to g^ u [3]: 

In this expression g{x^) is the determinant of —g^ v (x^) and g^ v is to be considered 
an external field rather than a dynamical variable. For any metric that is close to 
Minkowskian, i.e., g^ v = r]^ + h^ v and \J g(x^) ~ 1, one can easily see that the 
action differs from that of a free particle by | / d 4 iT M " '(x^h^x^) . Therefore, the 
interaction Lagrangian density is given by 

Cj = l -K u T^. (5.2) 

The Euler-Lagrange equations following from a Lagrangian with this interaction 
term are the same as implied by the equation of geodesic deviation only if is 
expressed in the LIF and T^ v is the classical stress-energy tensor of the system. The 
reason is that the Lagrangian ()5.2j) is not invariant under infinitesimal coordinate 
(gauge) transformations x'^ = x^ + ^, even though the action from which it is derived 
is a scalar quantity. Of course, this doesn't mean that the computed motions of 
particles depend on a particular gauge. The T^ u in Eq. (|5.2|) includes the effects of 
the classical, non-gravitational forces acting on the particle. In non-LIF gauges, e.g., 
the TT gauge, these forces can be significantly modified by the gravitational wave 
and the modifications must be taken into account in calculating particle motions. 

As a specific example, consider a system of particles held together by electro- 
magnetic forces. In principle, the action must therefore include the electromagnetic 
stress-energy, which is coupled to the gravitational field. In the LIF gauge gravita- 
tionally induced electromagnetic forces are smaller than the tidal forces and can be 
neglected. On the other hand, in a non-LIF, the solution to the Einstein-Maxwell 
equations include gravitationally induced electromagnetic forces that are compara- 
ble to the tidal effects and, therefore, they must be taken into account. (A simple 
example of such a system was treated by Boughn |S].) 

Dyson's treatment of elastic systems in the presence of gravitational waves JTj is 
based on the interaction Lagrangian ()5.2|) as expressed in the TT gauge, which does 
not constitute a LIF. Gravitationally induced modifications to the internal stresses 
must therefore be taken into account. Dyson does not do this, however, and concludes 
that "the response [of an elastic solid] depends on irregularities in the shear-wave 
modulus, and is strongest at free surfaces." As a consequence, Dyson's analysis 
implies that a self-gravitating, compressible perfect fluid (a system without shear), 
should not interact with a gravitational wave. But in fact the sun is a reasonable 
approximation to such a system at periods comparable to sound travel times across 
the sun and yet does couple to incident gravitational waves |T21 EE] • 

Similarly, in computing their transition rate in hydrogen, the Problem Book au- 
thors [TU] worked in the TT gauge but did not include electromagnetic stresses in 
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the interaction Lagrangian. Recall that in special relativity, the Minkowski metric 
can be written dr 2 = —7] liu dx tl dx l/ = dt 2 (l — v 2 ). This leads to the free-particle 
Lagrangian L = —m\J\ — v 2 . For a charged particle interacting with both an elec- 
tromagnetic and a gravitational field we replace r]^ by g^ v = i]^ v + h^ v and add the 
usual electromagnetic interaction term, giving 



-m 



yl — f]jkV^v k — hjkvH k — g$ + qv J Aj, (5.3) 



where $ and A are the electromagnetic scalar and vector potentials and we assume 
the h^ = condition of the TT gauge. For non-relativistic systems Eq. ()5.3|) can be 
expanded to lowest order in u J , 

L = —m + ^m(r]j k v j v k + h jk v j v k ) — g$ + qv j Aj. (5.4) 

Note that the hj k term in this expression is just the interaction Lagrangian of Eq. 
(|5.2|) . The Hamiltonian is formed through the standard prescription H = 7T a q a — L, 
where n a = dL/dq a and q a are the generalized momenta and velocities. Working to 
first order in hjk the canonical momenta are 



m(r] jk v k + h jk v k ) + qAj. 



To first order, the inverse of (f]j k + hj k ) is (f]^ — h^) where h? k = hj k . Then solving 
for v i gives 

v* = -(rj jk -h jk )(7r k -qA k ) 
m 

and the Hamiltonian becomes 

2m 

The constant m is not physically relevant and, as one sees, the Hamiltonian ()5.5j) 
follows the usual minimal substitution rule, tt — > 7r — qA. To implement Eq. ()5.5|) . 
TCj is identified with the quantum mechanical operator —ihd/dx^ and hj k is identified 
with the graviton raising and lower operators according to Eq. 14.21 For nonrelativistic 
systems, magnetic fields are much smaller than electric fields and so, in the Coulomb 
gauge, the vector potential A can be ignored. In this case the interaction part of the 
Hamiltonian operator becomes 

H - ^ h- 8 8 (5 6) 

2m 3 dxj dx k 

When this expression is substituted into the matrix element of Eq. ()4.3j) the result 
is a transition rate 16/25 times that of Eq. ()4.6j) and is the rate presented in the 
Problem Book, modulo numerical and normalization errors. 
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This factor can be accounted for by including the purely electromagnetic part of 
the Lagrangian density 

Cem = = --LgrtpFopFp, (5.7) 

I07T l07T 



where F^ = A v ^ — A^ v is the usual electromagnetic field tensor [14J. To lowest order 
in h? k this becomes 

C EM = -z^-V^F^F^ + ^-h^rTF^Fok. (5.8) 

The purely electromagnetic part of the Hamiltonian density then has the form 

Hem = (dC EM /8A jfi )A jfi - C EM - (5.9) 

It is straightforward to show that to first order in hP the total Hamiltonian becomes 
(see, for example, Chapter 57 of Schiff |3]) 

H = ^-h^-qA^-qA,) + J_ I 2£x + hj%Ek)d s x (5 1Q) 
2m o7r J 

where we have again invoked the Coulomb gauge and ignored terms involving the 
magnetic field so the electric field is given by Ej = A 0) j. Finally, the interaction part 
of the Hamiltonian can be written as 



h 2 8 8 1 r 
R I = — h jk —— + — / h jk E j E k d 3 x (5.11) 
2m 3 8x,8x k 8vri 3 K J 



h 1 , 8 8 

Cj 8x k 



since h? = hj k and to lowest order in hj k there is no distinction between upper and 
lower indices in the other fields. 

In this expression E^ is the sum of electric fields of the proton and the electron, 
E 3 p + E{ , so the integrand of Eq. (I5~TT) is (£» + E{) (E* + E k e ) . The integrals of E 3 p E k 
and E{E k e over all space are constants, independent of the location of the particles. 
Therefore, the contribution of these terms to the matrix element (see Eq. 14. 3|) between 
the Is and 3d0 states vanishes. The only electromagnetic terms in the interaction 
Hamiltonian that contribute to the transition rate are 

H EM ,i = ^h jk J (Ej,E k + E k Ei)d 3 x (5.12) 

where we again assume that hj k (r) is approximately constant in the region of signifi- 
cant electric fields.. Consider the first term 
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Here e is the electron charge, r is the field point, r e is the electron coordinate, and 
the proton is assumed to be located at r = 0. Via the addition theorem for spherical 
harmonics, |r — r e | _1 can be expressed as [T4"] . 

1 00 1 1 r 



47r ZJ ZJ -- - rhr YL(ff,&)Y f 

£=0 m=-l 



^ ^ 2£ + lr m lmK ,V} im{ 



Noting that the matrix element between the Is and 3d0 states will pick out only the 
I = 2, m = term, we are thus left with 

1 Airr 2 

V —> = -^Y* {9, 0)F 2O (^e, <f>e), (5.14) 

where r< (r>) is the smaller (larger) of r and r e . Substituting Eq. (j5.14j) into Eq. 
(I5.13jl and integrating over all space yields, after a tedious calculation, 

2to „„j , 



/ d?xE;Et = J SxElEl = -\j SxElK = 1} - (5.15) 

The integrals vanish for j ^ k and, therefore, the two terms in Eq. (|5.12|) are 
equal. Substituting Eq. (|5.15jl into Eq. (|5.12j) and adding the result to the electron 
Hamlitonian (see Eq. 15.6)1 yields the interaction part of the Hamiltonian. Another 
tedious computation via Eq. (|4.3J) gives, finally, a transition rate of 

T = P '^= 5x2 i3ft c 5 =-36M*-' (5 - 16) 

which is precisely the same as the LIF field theoretic and semi-classical results (|4.6jl 
and (1231) . 

From this analysis we conclude that one can work with an interaction Hamiltonian 
in any gauge as long as all the revelent interactions are included in the stress-energy 
tensor. However, it is also clear that in this case working in the LIF leads to a much 
simpler analysis. For classical systems, it is undoubtedly also possible to work in 
non-LIF gauges; however, in the case of elastic media considered by Dyson [7j, for 
example, the best way to do this is not immediately apparent. 



6 Graviton Ionization cross section of Hydrogen 

In RB, an important step was to calculate the ionization cross section for hydrogen. As 
we saw in §21 the graviton absorption cross section was o ~ £ 2 pl . That all dimensions 
other than the Planck length drop out of the cross section is at first surprising, but 
can be understood simply as follows: The classical cross section for a system of mass 
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and frequency u is roughly 0||2| cr ~ Gm£ 2 uj/c 3 . Assuming the Nicholson- 
Bohr quantization condition, the angular momentum for such a system near its ground 
state is L ~ ml 2 uj ~ ft. Thus 

a ~ ^ = 4 Z « 10- 66 cm 2 , 
cr 

and we see that the Planck-length-squared cross section is solely a result angular 
momentum quantization. 

Because in RB we decided on the gravitational analogy of the photoelectric effect 
as a method for detecting gravitons, it was necessary to compute the gravito-ionization 
cross section for hydrogen in the ground state. One expects it to be similar in mag- 
nitude to the above; however the ordinary photoionization cross section does have 
large numerical factors "of order unity" and is also strongly dependent on photon 
energy. Therefore, in this section we compute the ionization cross section for gravi- 
tons with energies in the range 13.6eV <C E <C 2.5 x 10 4 eV, energies for which the 
non-relativistic Born approximation is appropriate. In fact, as detailed in RB, there 
are many strong astrophysical sources of gravitons in this energy range. 

As in JJJ we take the interaction Hamiltonian in a LIF to be 

H = -m e u} 2 hxix k ejkex i pi(kix l — ut) + c.c. (6.1) 

and compute the matrix element between an initial hydrogen ground state, and 
a plane wave final state (the Born approximation), i.e., 

where the plane wave is box normalized with dimension L. The transition probability 
per unit time between these two states is given by the golden rule, 

r = ^p(k)\(f\H\i)\* (6.3) 

where k is the wave number of the emerging electron and p(k) is the energy density 
of final states: 



Pik) = (6.4) 



As in JJl we average \(f\H\i)\ 2 over all directions using Eq. (|3.12j) . For an incident 
GW of amplitude h this result is 



\(f\H\i)\ 



2 3 • 2 n 7r h 2 uj A m 2 e a 7 (a A k A ) 



L^l + a 2 k 2 Y 



Gravitons 



17 



Eqs. fTHD and (E3D then give 

_ 3 ■ 2 10 h 2 ^m\a xl k h 
~ 5 h 3 {l + a 2 k 2 ) s ' 

By definition, the ionization cross section is a = Thu jT where T = h 2 is the 
GW flux from m Thus, 



_ 3 • 2 13 tt Gu 3 m 3 a n k 5 
a ~ 5 c 3fr 2 (l + a 2 P) 8 - (6 ' 5) 

uj can be eliminated from this expression by using the Einstein photoelectric relation 
(conservation of energy), which requires that the incident graviton energy equal the 
sum of the electron binding energy and the kinetic energy of the emerging electron, 
or huj = e 2 /2a + k 2 /2m e . Hence, 

_ 3 ■ 2 10 7r (kaf Gh 
a 5 (1 + k 2 a 2 f c 3 ' 1 ' ' 

For ak » 1 the dependence of this result on final electron momentum is in fact 
identical to that of the ordinary photoionization cross section. We also see that the 
ionization cross section is, modulo a dimensionless factor, equal to the Planck length 
squared; for energetic gravitions with Tiuo » 13.6eV, however, the dimensionless 
factor can be quite small. Smolin[ro] used a field theoretic argument to estimate 
the ionization cross section of a bound electron and obtained a result that is also 
proportional to £ 2 Pl . However, the dependence on ka is quite different, with the cross 
section increasing rapidly for large ka rather than decreasing rapidly as indicated in 
Eq. (IFTTIjl . 

The differential ionization cross section for linearly polarized gravitons is also not 
difficult to compute and we state it here for completeness. For one polarization we 
obtain 

^ n = 32 ' 25 (TTI& ? sin22 ^ in49 ? i; 

letting; sin 20 -> cos 2 20 gives the other polarization. 



7 Conclusion 

As discussed in the Introduction, the original motivation for this paper was Dyson's 
conjecture that a single graviton could never be detected in the real universe. In RB 
we employed the above gravito-electric cross section to show that if one is limited 
only by the mass, energy content, and age of the universe, one can design a highly 
idealized gedanken experiment that could detect some gravitons in the lifetime of the 
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universe. As soon as one begins to consider detector physics and background noise, 
though, detecting even a single graviton becomes impossible. In that sense, Dyson's 
conjecture appears correct. 

Although one might argue that the detailed calculations presented here are not 
entirely necessary to address Dyson's conjecture, it has become clear that the physics 
of gravitational wave-matter interaction, in particular the gauge properties of the 
interaction Hamiltonian, present enough subtleties to catch even experienced practi- 
tioners off guard. As this will undoubtedly happen from time to time in the future, 
we feel it is of some importance to elucidate this matter as clearly as possible. 
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